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1 Introduction 

The bipolar expansion of the inverse distance between 
two points and its applications in the study of chemi¬ 
cal bonding has been widely studied in the literature 
by different methods [lj|^, as being very suitable for 
approximating the interaction potential between vari¬ 
ous kinds of molecules in terms of their multipoles. It 
is usually used in a four-region form, which imply the 
decomposition of the integration space in regions with 
complicated boundary-surfaces: a nonoverlapping re¬ 
gion and three overlapping regions. As many modern 
numerical calculations used in quantum chemistry are 
based on multipolar approximations [^[^, the con¬ 
vergence of the bipolar expansion is very important. 
The recent conclusions of the numerical studies from 
Refs. are not encouraging at all. They reveal 

that the convergence in the overlapping region is con¬ 
ditional and that the Laplace equation is not satisfied 
either termwise or pointwise. Nevertheless, the au¬ 
thor of hopes in ’’the possibility of clarification 
through generalized functions”. 

In this paper, we are going to analyze new forms 
of the bipolar expansion, deduced from the multipolar 
formalism of Refs. m and calculate the overlap cor¬ 


rection to the long-range electrostatic energy for two 
simple systems; a hydrogen atom and a proton and 
two hydrogen atoms in their ground states. Using this 
formalism we obtain well-known results 
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and this 

encourages us to assert that our expansions are ” good 
enough” for studying more complicated systems. In 
contrast to other papers which are studying the bipo¬ 
lar expansions, we start from the multipole series of 
the electrostatic interaction energy and compare it 
with its definition as can be found in textbooks. Our 
results contain well-behaved functions as well as gen¬ 
eralized functions, therefore their convergence can be 
discussed only in the weak or distributional sense. In 
our approach, using the results of Ref. , we are able 
to identify the generalized-functions terms with the 
contribution of the mean square radii and to prove 
that they describe the overlap (or contact) interac¬ 
tion. 

The paper is organized as follows: in the next section 
we justify four equivalent forms for the bipolar ex¬ 
pansion of the inverse distance, valid in all regions of 
space. Next, we use the perturbation theory to cal¬ 
culate the first-order correction to the electrostatic 
interaction energy for two simple physical systems. 
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neglecting the exchange. The end section is devoted 
to conclusions. 

2 New Forms of the Bipolar Expansion 

We consider two charge distributions, one of which is 
centered at the origin and is described by the vector b 


and the other centered at a distance R from the origin 
and described by the vector a, as shown in Fig. 1. We 
are concerned with the inverse distance between the 
points A and B. We shall justify here the following 
four equivalent equations: 


z 





Figure 1: Coordinates for the bipolar expansion 
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where we have used the following notations: 

• Vi is the i-th component of the vector r and 
r = |r|, , r = r/r 

- Yirai^) is a spherical harmonic and Yirair) = 
f^Yiraii) is a solid spherical harmonic ( |l2[ ) 

- dri = ^ means partial derivative with respect 
to the i-th component of the vector r 

- V is the gradient operator and A = V • V is the 
laplacian 


- aj? = dr_^, A" = Ar ■ ■ ■ Ar 

n times n times 


The summation over dummy indices is understood. In 
the following, for shortness, we shall habitually drop 
the limits in summations. 

According to Ref. , the two charge distributions 
are described by the expansions: 


(/+/', is a Gaunt coefficient 13 
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where p(r) is the charge density, (5(r) is the Dirac 
delta function, Sim{r) = 

two sets of mean square radii ^ are calculated 

with respect to the two centers of the charge distri¬ 
butions: 

= Vsfi / 

The electrostatic potential produced by the sec¬ 


ond charge distribution is: 

A ('^'1 — + 1) / ^TT 

(P 2 { ) - 2^n!(2/ + 2n + 1)!! V + T 

l,m,n 

( 8 ) 

After integrating by parts, one obtains the follow¬ 
ing expression for the electrostatic interaction energy 

IT = y (frpi{r)(j)2ir) : (9) 
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the mean square radii 0 and interchange sums and 


( 10 ) 


W = f dh /da pi(b) P2(a) ^ 

1,1' ,m,m' ,n,n' 




2n+n'n\n'\(2l + 2n + 1)!!(2^' 2n' + 1)!! 


y;:„(b)y;^,(a)A-+-V,^(-VR)y,,^,(-VR)y. 


( 11 ) 


If we compare the above equation with: 


W = 


db 


da 


Pi(b) p 2 (a) 

lb — a — R| 


we easily obtain the expansion ([^ . Eq. Q can be eas¬ 
ily verihed in the Fourier space: if we take the Fourier 
^]^2) transform f dR and make the summations over 
the indices n and n' we obtain: 


e*k(b-a) = jK^?^)jT'(fca)>):n(b)k)r^.(a)yz™(k)yzw(k), 

1,1' ,m,m' 


( 13 ) 


4 



















where 


ji{x) - Am ^ n\2^(2l + 2n + l)\\ 

are the spherical Bessel functions. Thus, we have 
obtained the well-known Rayleigh expansion of the 
plane wave. Eq.Q is a slightly modified form of 
Eq.Q, obtained by using the addition theorem of 
the irregular solid spherical harmonics from Eq. 
(7.4), which we transcribe here in a slightly modified 
form: 


After we use the dehnition of the Pochhammer sym¬ 
bol: {a)n = > where r(a:) is the Euler Gamma 

function, the above equation becomes: 

VR,)Fz/m/(—V r,)— = 

{l + l' m + ( 15 ) 


and our Eq.(2) follows immediately. 

Now, in Eq.(|10[) we pass to the Cartesian tensor 
formalism, using Eq.(16) of Ref. 10 , which can be 


written in operatorial form as follows: 
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We obtain the following expression for the electro¬ 
static interaction energy: 


where 


W = 


E 


(-l)'(2i + l)!!(2f+l)!! ^(1) ^(2) 

, 2^+’^'n!n'!/!r!(2Z + 2n + l)!!(2f' + 2n' + l)!! -h -T' 


,n,n 

/\n+n'i^, Q. D. D. 

Ojj . . . . . . Uj^, 


1 

R’ 


of Ref. [^), as follows: 
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are the Cartesian components of the mean square ra¬ 
dius of order n and multipolarity I, and the summa¬ 
tion over dummy indices is understood. If we intro¬ 
duce in Eq.(|17[) the above expression for the Carte¬ 
sian mean square radius and compare with Eq.(12), 
we obtain Eq. (§. Eq.Q can be obtained by writing 
the contractions of the totally symmetric tensors in 
terms of the compressed forms of the tensors (Eq.(2) 
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3 Calculation of the Long-range Electrostatic 
Energy 

3.1 The interaction between two hydrogen 
atoms 

In this section we consider the long-range electro¬ 
static interaction between two hydrogen atoms, in 
their ground states. We shall neglect the exchange 
and calculate the first order correction to the energy 
by using the non-degenerate perturbation theory. As 
the charge density of a hydrogen atom centered at the 
origin is: 

p{r) = e6{r) — e6{r — a), 


where e > 0 is the charge of the proton, the mean 
square radii of the two interacting systems (calculated 
in their own reference systems) are: 



^2n' 

'Vm' 




Introducing the above expressions in Eq.(lO), we ob¬ 
tain the following expression of the perturbation op¬ 
erator: 
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Note that, if we place the x-axis along R and explic¬ 
itly write the first terms of the above expansion, we 
obtain the r.h.s. of Eq.(3) from |14| , supplemented 
by an infinite number of point-like terms (which con¬ 
tain Dirac delta functions and their derivatives). This 
point-like terms are the contribution to the energy op¬ 
erator of the mean square radii: for each multipolarity 
order, there exists an infinite number of mean square 
radii, which contribute only to the overlap (contact) 
energy. Considering the contribution of the mean 
square radii is equivalent to considering the correct 
form of the bipolar expansion in the interpenetrating 
region. 

As pointed out in Ref. [^, the first consequence 


of using the correct form of the bipolar expansion in 
the overlap region is that the first order correction 
to the energy no longer vanishes (in Ref. 14 , where 


the contribution of the mean square radii have not 
been considered, the first-order correction to the in¬ 
teraction energy between two hydrogen atoms in their 
ground states is zero). 

The calculation of the first-order correction to 
the energy can be done either in the coordinate 
space (using the operator ®) or in the recipro¬ 
cal (Fourier) space. We shall use here the second 
method. We take the Fourier transform W (a, b; k) = 


/dR 


..zkR 


iy(a,b;R) and, after summing over the 


indices n, n', we obtain: 
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Due to the spherical symmetry of the ground- give contribution to the first-order correction to the 
state wave function of the hydrogen atom, only the energy: 
monopole mean square radii [I = I' = m = m' = Q) 
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is the hydrogen ground-state wave function, ao is the 
Bohr radius. All the integrations can be easily done 
and, after taking the inverse Fourier transform, one 
obtains: 


eC)(r) = ;^ exp(- —). 
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where Eq = —e^/(2oo) is the ground-state energy 
of the hydrogen atom. This is the result obtained 
in Eq.(31) of Ref. 11 , by the direct integration of 
the Coulomb potential, when the exchange is ne¬ 
glected. As in our calculation this term have been 
obtained from the contribution of the monopole mean 
square radii, we conclude that it has nothing to do 
with the polarization of the atoms. The effect of the 
polarization (contributions from the dipole operator, 
quadrupole operator etc.) appears in the second order 
of the perturbation theory. 


3.2 The interaction between a hydrogen atom 
and a proton 

The unperturbed system consists of a hydrogen 
atom in its fundamental state and a proton. We 
study their electrostatic interaction, neglecting the 
exchange (that is, we suppose that the electron does 
not jump from one proton to the other). The charge 
mean square radii of the two interacting systems are: 
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The perturbation operator is: 
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and its Fourier transform: 
tF(a;k) = 

l,m 

dvre^ 

One obtains: 


(25) 


^(b(k) = J da^p^°>{a)W{a■k)^|;^^\a.) = 

dvre^ / 16 \ 
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and after taking the inverse Fourier transform we get 
the first-order correction to the energy: 


E^^'){-R) = ^e + 1^ (27) 

R Voo J 

This is the well-known result for the first-order cor¬ 
rection to the energy, when the exchange is neglected 
(Eq. (19b) of Ref. [^). As in the preceeding sub¬ 
section, this correction to the energy comes from the 
contributions of the monopole mean square radii of 
the hydrogen atom, so it has nothing to do with the 
polarization. 


4 Discussion and Conclusions 


Our main results 0 - ( |10[ ) contain 

series of the type Ez,m,n[-• • ]^rn(-V)A";^ = 

-4vrX)z,m,n[-• Such series of 

Dirac delta functions were studied both by math¬ 
ematicians and physicists, in connection with the 
asymptotic analysis, generalized solutions of differen¬ 
tial equations, multipole analysis, quantized field the¬ 
ories, etc. |16| - [23| . Starting from the need of multipole 
series in different branches of physics, the mathemati¬ 
cians extended the Schwartz’s space of distributions 
to different spaces of ultradistributions (polynomial, 
tempered, with exponential growth, of compact sup¬ 
port, Hermitean, etc), in which some series of Dirac 


delta functions are convergent. Without entering into 
the mathematical details of these theories, we wish 
to underline here some important aspects regarding 
these distributional series: 

• our formalism make sense only for finite or 
rapidly descending charge densities, so that the 
mean square radii Eq.0 should be finite. 

• the distributional series which appear in our re¬ 
sults are convergent at least as polynomial ultra¬ 
distributions. Finding the largest space in which 
our results are convergent remains an open prob¬ 
lem. 

• these infinite series of Dirac delta functions are 
not localized (supported) at the origin ( [Isp^ ). 
If we regard them as functionals applied to some 
analytic function cj), we obtain: 


OO OO 

(28) 

n=0 n=0 


so they depend not only on the values of 4> at 0, 
but on the values of cf) in a certain neighborhood 
of 0. In our multipole formalism these series are 
precisely connected to the finite or inhnite over¬ 
lap of the charge densities. 

• When we calculate the interaction energy by us¬ 
ing Eq. 0, if the charge density /9i(r) is an an¬ 
alytic function, we obtain a convergent power 
series. Let us consider, for example, the in¬ 
teraction between a uniformly charged sphere 
/?i(r) = Ci6{a — r) and a gaussian charge den¬ 
sity / 02 (r) = . From Eqs.0, 0, one 

obtains: 

6a 15 

(29) 

which are the first two terms of the Taylor series 








of the exact result: 
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where 7 is the incomplete Gamma function. 


• all the derivatives which appear in our formu¬ 
las should be taken in distributional sense. Let 
us consider the interaction between a uniformly 
charged sphere /Oi(r) = Ci6[a — r) and an expo¬ 
nential charge density P 2 {^) = If we 

use Eq. Q for the electrostatic potential of the 
sphere and calculate W = f dr (/>i(r)p 2 (r — R), 
when we integrate by parts we should apply the 
operator A” to the exponential function, with¬ 
out neglecting the Delta-type terms: 
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Note that the series which contain a2"'+3 from 
the above equation (lines 2, 4 and 6 ) should not 
be viewed as power series in ”a” (because their 
coefficients are generalized functions) and are 
not useful in numerical calculations. For ii = 0 
one obtains: 


— Ar 
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For n = 0, 00 , we obtain a series of delta func¬ 
tions and their derivatives, which should be 
summed. We can do this by using the Fourier 
transform method: 
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which coincide with the result obtained by di¬ 
rect calculation. 


Using the multipolar formalism, as presented in 
Refs. [8,10, we have written four compact equiva¬ 
lent forms for the bipolar expansion of the inverse 
distance, which are valid in all the regions of config¬ 
uration space. In the nonoverlapping region, these 
expansions coincide with those already known (Hi. 
The equivalence of our results in the interpenetrat¬ 
ing region with the other results from the literature 
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remains an open problem. Using the first order per¬ 
turbation theory, we have calculated the overlap cor¬ 
rections to the long-range electrostatic interaction en¬ 
ergy for two simple systems, neglecting the exchange. 
We have shown that, considering the contributions 
of the mean square radii to the interaction energy 
is equivalent to considering the correct form of the 
bipolar expansion in the interpenetrating region. 
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